Superflow-Stabilized Nonlinear NMR in Rotating ^He-B 
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Nonlinear spin precession has been observed in "^He-B in large counterflow of the normal and 
superfluid fractions. The new precessing state is stabilized at high rf excitation level and displays 
frequency-locked precession over a large range of frequency shifts, with the magnetization at its 
equilibrium value. Comparison to analytical and numerical calculation indicates that in this state 
the orbital angular momentum L of the Cooper pairs is oriented transverse to the external magnetic 
field in a "non-Leggett" configuration with broken spin-orbit coupling. The resonance shift depends 
on the tipping angle 6 of the magnetization as to — ujl = (f2s/2ti;L) (cos 6 — 1/5). The phase diagram 
of the precessing modes with arbitrary orientation of L is constructed. 

PACS numbers; 67.57.Fg, 47.32.-y 
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The study of nonlinear NMR response in superfluid 
■^He started with the discovery of the Brinkman-Smith 
(BS) mode in ^He-A § and ^He-B [|| phases. These be- 
came the classic examples of nonlinear spin resonance in 
magnetically ordered superfluids. Their observation by 
OsherofF and Corruccini |^ opened the road to the dis- 
covery of new resonance states in '^He-B, such as space- 
coherent precession within a Homogeneously Precessing 
Domain (HPD) Q or the newly found stable modes 
where the magnitude of the precessing magnetization dif- 
fers from its equilibrium value Q. An example of the 
latter is the family of half-magnctization modes (HM), 
where the processing magnetization equals one half of 
the equilibrium value All of these are stable dy- 
namic order parameter states and nonlinear solutions of 
the Leggett-Takagi spin dynamic equations. In principle, 
such states are similar to Q-balls, the coherent soliton- 
like nontopological states of relativistic quantum field 
theories, whose frequency and stability are determined 
by the conservation of the global charge, say, the bary- 
onic charge Q. In •^He-B spin dynamics the role of the 
global charge is played by the projection of the spin Sz in 
the direction of the magnetic field H , which determines, 
in part, the NMR frequency shift. 

In these resonance modes the distinguishing factor is 
the orientation of the orbital angular momentum L of 
the Cooper pairs. In the BS and HPD modes, L is ori- 
ented along the applied magnetic field H via the spin- 
orbit (dipole) coupling, i.e. L \\ H . In contrast, the HM 
modes form with L oriented spontaneously perpendicu- 
lar to H . So far, in the nonlinear regime the orientation 
of L has not been controlled by external means. Here we 
use vortex-free counterflow of the normal and superfluid 
components in a rotating container, to orient L along the 
flow direction. At high rotation velocity 17, the orienting 
effect on L from the flow far exceeds that from the dipole 
coupling. If the external magnetic field is oriented along 



the rotation axis {H \\ fi), one can then study the un- 
usual situation, when L H. By measuring the tipping 
angle of the precessing magnetization M as a function 
of the applied frequency shift, we identify a new nonlin- 
ear resonance mode which greatly differs from the classic 
case of L \\ H . In the linear regime the condition L L H 
has been realized in earlier measurements in the parallel- 
plate geometry or in the presence of counterflow 

Experiment. — Our cw NMR setup in the rotating nu- 
clear demagnetization cryostat has been described in 
Refs. [^IJq) . The "^He NMR sample is contained in a quartz 
glass cylinder with a radius R — 2.5 mm. The measure- 
ments are performed at fixed frequency Wrf/27r = 688 
kHz, using a linear field sweep centered around a Lar- 
mor field value of H = 21.2 mT, with a homogeneity 
AH/H = 2 • 10~* over the sample volume. The signal is 
read with a lock-in amplifier, such that the component 
in phase with the excitation field Hj-f is called dispersion 
(oc Mx = M± cos (/)) and the out-of-phasc component ab- 
sorption (oc My = M± sint/)). The measuring range com- 
prises counterflow velocities flR < 6.5 mm/s, rf fields 
up to 0.03 Oe, temperatures (0.7 — l)Tc, and pressures 
0—12 bar. 

At low excitation amplitude (i/rf ^ 0.003 Oe), when 
the tipping angle is a few degrees, the NMR response 
of ^He-B exhibits linear behaviour: The line shapes of 
the absorption and dispersion signals are independent of 
sweep direction, and the signal amplitudes increase lin- 
early with excitation level. In vortex-free counterflow at 
sufficiently high velocity {n > 0.5 rad/s) the NMR ab- 
sorption maximum is shifted from the Larmor value Jp|, 
as shown by the NMR spectrum in the inset of Fig. ^ 

At high excitation levels the absorption and dispersion 
signal amplitudes increase faster than the rf field, and be- 
come highly asymmetric (Fig. The nonlinear behavior 
becomes most pronounced while scanning the field in the 
upward direction towards the Larmor value until finally 
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FIG. 1. NMR response while the apphed magnetic field 
H is swept at constant rate H during continuous rf irradia- 
tion at fixed frequency lo^{ {Q, = 2.5 rad/s, P = Obar, and 
T ~ 0.88 Tc). Up sweeps are marked with thicker lines and 
down sweeps with thinner lines. Solid lines correspond to 
i/rf = 0.02 Oe and H = 0.017 Oe/s (Exp. 1), dashed lines 
to Hrt = 0.027 Oe and H = 0.034 Oe/s (Exp. 2). Note that 
at higher excitation lower absorption is measured for compen- 
sating the relaxation. (Insets) Line shapes of the signal in the 
linear regime at low exitation level (SI = 1.3 rad/s, P = 2.0 
bar, T = 0.90 Tc, and H,{ = 0.003 Oe). 

an abrupt jump appears from the precession at large tip- 
ping angle 9 to the linear NMR regime with small 9. If 
the excitation amplitude is increased, the jump usually 
moves to higher field. The field sweep in the opposite di- 
rection has significantly different shape and the regime of 
large tipping angles is not entered. The magnitude of the 
counterflow plays a crucial role (Fig. ||): With decreasing 
n both the maximum tipping angle and the range of fre- 
quency shifts quickly decrease. This means that the new 
state appears only at high counterflow velocities above 
the textural transition in which L is deflected into the 
plane transverse to f2 in a significant part of the cross 
section of the sample cylinder . 

We interprete these observations in the following man- 
ner. For fixed orientations of S and L, the resonance 
frequency shift is determined by the dipole interaction. 
During the field sweep these vectors defiect from there 
equilibrium positions and the frequency of the resonance 
absorption changes. If the sweep is performed in a suit- 
able direction, it may become possible to create a state 
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FIG. 2. Transverse magnetization component Mj_ as a 
function of field H, at different counterflow velocities in vor- 
tex-free rotation (listed in terms of Q and the dimensionless 
variable u, calculated at the outer perimeter of the sample 
cylinder from Eq. (|^). Only upward sweeps of the field H 
are shown, down sweeps depend less on Q and are similar to 
those in Fig. ^. The temperature increases slowly during the 
measurement from 0.88 Tc at the highest velocity to 0.92 Tc 
at the lowest velocity. The normal state magnetization Xn^^ 
corresponds approximately to 200 units on the vertical scale 
(P = Obar, H,{ = 0.027 Oe, and H = 0.017 Oe/s). 

where the resonance frequency stays locked to the ex- 
ternal excitation frequency. In this case the increasing 
deflection of S causes the transverse magnetization Mj_ 
to increase, which is observed in the experiment during 
a sweep towards the Larmor value as an increase in the 
dispersion and absorption signals. The tipping angle in- 
creases continuously with the fleld sweep as long as the rf 
pumping is sufficient to compensate for relaxation, which 
increases as the deviation of S from its equilibrium ori- 
entation increases. Finally the mode collapses, in a first 
order transition between two different dynamic order pa- 
rameter states. This behaviour resembles that of an an- 
harmonic oscillator in forced oscillation. The remarkable 
feature of ^He-B is that the rigidity from the order pa- 
rameter coherence makes the superfluid to behave like a 
single oscillator. 

Classification of modes. — An analytic description of 
spin precession with arbitrary orientations of S and L can 
be constructed, if we neglect magnetic relaxation and the 
interaction with the excitation field. The orientation of 
the orbital momentum, below denoted by the unit vector 
I — —L/\L\, is fixed by the balance between its interac- 
tions with the counterflow and with the precessing spins 
via the dipole coupling. The former is written as 

Fcf = -^PaiiVs - Vn) ■ i)^ . (1) 

Here pa — Ps± ~ Ps\\ is the superfluid density anisotropy, 
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caused by the magnetic field jlj] . The relative magnitude 
of the counterflow and dipole energies is conveniently ex- 
pressed in terms of a dimensionless velocity: 



cos 9 



(2) 



where ftsiT, P) is the characteristic ■^He-B frequency and 
a measure of the dipole energy. Spin precession is simpli- 
fied in the high- field limit, when the dipole term is small 
compared to Zeeman energy, —S-jH, and can be consid- 
ered as a perturbation. In zero order perturbation theory 
one has precession with the Larmor frequency, lol = . 
A first order correction gives the frequency shift from 
the precessing frequency lo: lo — ujl — —dFu/dSz, where 
Fo is the dipole energy, averaged over the period of the 
precession. Here we consider only the case when the pre- 
cessing spin has its equilibrium magnitude S = xbH/j- 
For arbitrary orientation of the orbital momentum, the 
time-averaged dipole energy Fd can be written as 
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l^nliisj, - 1 -I- icos<i>(l + -I- + 
(1 - szfil - Uf + (1 - _ + cos$)] . (3) 



The notations are: = Szl\S\ = cos 6*, 1^ is the projec- 
tion of the orbital momentum I on the direction of the 
magnetic field if, and the angle $ is a soft variable re- 
lated to the ■^He-B order parameter. The energy Fu is 
also to be stationary with respect to $: dFujd^ = 0. 

Thus by varying the dipole energy with respect to the 
spin 5"^ (the analogue of the global charge) one obtains 
the precessing modes as a function of the frequency shift 
and of the other global charge L^, which is kept fixed 
because of orbital viscosity. Omitting the limiting cases 
6* = and = we get three modes of precession 
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2) cos$ = l, 

3) cos<I> = -l, 
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Here w is the dimensionless frequency shift: 
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In Fig. 1^ the Iz — Sz phase diagram is shown with the 
stable regions for each of the three modes. In the absence 
of counterflow, u = 0, the dipole coupling orients the 
orbital momentum I along the magnetic field H in modes 
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FIG. 3. (Left) Phase diagram of the precessing states for 
arbitrary orientations of L and S. (Right) Dependence of the 
frequency shift w on the tipping angle d for a few fixed values 
of L (marked for each line at its end point) . 



1 and 2 and opposite to the field in mode 3. Then mode 
1 becomes the BS state with zero frequency shift in the 
range — ^ < < 1 and Iz = I- Mode 2 reduces to the 
HPD state with —l<Sz<—\,lz = ^, and a frequency 
shift vs tipping angle dependence as w = — 2 — 8sz, while 
mode 3 transforms to the so-called HPD (2) which has not 
been seen experimentally JTst . In the generalized phase 
diagram of Fig. ^ with nonzero u, we retain these names 
for the regions in which their respective u — modes lie. 
In large counterflow linear NMR at small 9 and nonlinear 
NMR at large are located in Fig. ^ on the Iz — axis 
and belong to the same class of BS states. During an 
up sweep the precession moves continuously from small 
to large 0, until ultimately magnetic relaxation causes 
an instability and a first order transition takes 6 back 
into the linear regime. This is similar to a gas-liquid 
transition, where no symmetry break occurs. 

Limit of large counterflow. — In rapid rotation u be- 
comes large (u ^ 15) and the orbital momentum is rigidly 
forced into the transverse plane over most of the cross sec- 
tion of the sample cylinder: Iz = 0. The HPD state does 
not exist in this limit, while the BS and HPD(2) modes 
have the frequency shifts 
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In Eq. (g|) the linear regime of small 9 corresponds to 
Sz ~ ^ and w — i or uj — ujl = {2/5)fl'^/ujL. In the 
nonlinear regime at large 9 the frequency shift has a pos- 
itive slope as a function of the longitudinal magnetiza- 
tion: duj/d cos 9 > 0. This suggests that in free preces- 
sion in a pulsed NMR measurement homogeneous pre- 
cession should become unstable and break into domains. 
In contrast, in continuous rf excitation the phase of the 
precession is locked to that of i/if and no instability oc- 
curs. The HPD(2) shift in Eq. (m has a negative slope, 
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du/d COS 6 < 0, and, since the conventional HPD state is 
unstable in counterflow, HPD (2) is thus the only inher- 
ently stable mode in large counterflow, with spontaneous 
phase-coherence in free precession. However, the HPD(2) 
mode displays large Leggett-Takagi relaxation because of 
its large deviation from the Leggett configuration. Pre- 
sumably in the T — > limit, where relaxation vanishes, 
the HPD(2) mode might become observable. 

Comparison with experiment. — In Fig. ^ we plot the 
tipping angle 9 from Eq. (||) as a function of the fre- 
quency shift w, along with the two measured NMR re- 
sponses from Fig. |l|. The frequency shift, at which the 
new mode collapses in the experiment, is determined by 
the relaxation processes. At the pressure of 12 bar, re- 
laxation is less and the new state is often stable during 
the upward sweep until above the Larmor field value, 
i.e. to negative frequency shifts, in agreement with the 
right panel of Fig. ||. Overall, we regard the agreement 
in Fig. U as satisfactory, if we allow for two experimental 
difhculties. 

One uncertainty arises from assigning the proper value 
to temperature. Temperature is measured in the linear 
NMR regime at low rf level, using the fact that the coun- 
terflow absorption maximum is then centered at w = 3 
(inset of Fig. |l|). This frequency shift determines the 
Leggett frequency VLb(T,P) in Eq. ^ and, once cali- 
brated Q , can be used as a thermometer with a sensitiv- 
ity better than 0.001 T^. In the NMR response at high rf 
level the only feature, which qualifies for thermometry, 
is the maximum of absorption during the down sweep 
(Fig. |l]). Its location is not exactly w — 'd. In fact, our 
numerical simulations suggest that it is w < 3. This 
means that we ascribe a higher temperature and smaller 
r^B to our data than actually would be the case, which 
explains why the measurements are shifted to the right 
in Fig. 

Another factor is the heating of the sample by the ab- 
sorbed rf power. The sample cylinder is connected with 
a narrow channel to the refrigerator to prevent vor- 
tices from leaking into the NMR volume. The thermal 
resistance of the channel can lead to unaccounted tem- 
perature rise and distortion of the line shape during the 
field sweep. The heating is less important with a faster 
rate H of sweep. All data in this paper were measured at 
large rates so that the responses for up and down sweeps 
agree in their overlap region at small tipping angles, in- 
dicating that the temperature is approximately constant 
during the field cycling. The overall uncertainty in tem- 
perature we estimate to about AT = -1-0.02 Tc- This is 
sufficient to explain the difference between measurement 
and Eq. (|) in Fig. § 

Numerical simulation. — We have supplemented the an- 
alytic description with direct numerical solution of the 
Leggett-Takagi equations, by calculating the response of 
the spin-dynamic variables in the time domain for the 
spatially homogenous case during continuous rf excita- 




FIG. 4. Tipping angle ^ as a function of the normalized fre- 
quency shift w. Comparison of Eq. to the measurements 
of Fig. |l|. A numerical solution is also shown with param- 
eter values as in the experiment, except for H^f = 0.014 Oe, 
T — 0.95 Tc, H = 0.02 Oe/s, and spatially homogeneous coun- 
terflow at 0.70 cm/s. Upward sweeps of the field are plotted 
with thick and downward sweeps with thin lines. 

tion. Counterflow gives rise to the orientational energy 
Eq. (|l|) and to an additional torque, which contributes 
to Leggett-Takagi relaxation: T = 5Fcf/59, where 69 
is an infinitesimal 3-dimensional rotation in spin space. 
The experimental parameters (T, P, H, H, H^i, fl) are ad- 
justed to match the experimental conditions. A typical 
result is shown in Fig. |^. It agrees surprisingly closely 
with Eq. (|^), demonstrating that the effects from relax- 
ation and rf irradiation to the frequency shift are small. 
It also reproduces the shape of the measured NMR re- 
sponse, showing that the main difficulty in the compari- 
son is the shift of the measured data to a higher temper- 
ature. In general, the simulation result is found to move 
closer to Eq. (||), when T, Hif, or 17 are increased. These 
are all changes, which help to boost either the value of u 
or improve the compensation for relaxation, and thus en- 
hance the stability of the new mode towards larger fields 
during the field sweep. 

In conclusion, we have observed a highly nonlinear 
NMR response of "^He-B when the direction of the or- 
bital momentum is fixed by large counterflow. This new 
state of precession has been identified as a Brinkman- 
Smith mode in a general classification scheme of states 
with fixed direction of L. Other branches of the phase 
diagram can perhaps be found in experiments where the 
direction of L is fixed with solid walls or by tilting the 
magnetic field towards the flow direction. 
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